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Executive Summary 

Gamma-ray bursts (GRBs) are the most brilliant events in the universe. The intrinsic lumi- 
nosities of the bursts span more than five decades. At first glance, therefore, these events would 
hardly seem to be promising standard candles for cosmology. However, very recently, relations 
between the peak energy .Epeak of the burst spectrum in vF^, the isotropic-equivalent energy Eiso 
of the burst, and the radiated energy E-y of the burst - all in the rest frame of the burst source - 
have been found. In a way that is exactly analogous to the way in which the relation between the 
peak luminosity and the rate of decline of the light curve of Type la supcrnovae can be used to 
make Type la supernovae excellent standard candles for cosmology, so too, the relations between 
£^peak) Eiso, and point toward a methodology for using GRBs as excellent standard candles for 
cosmology. In addition, GRBs occur over the broad redshift range from z = 0.1 to at least z = 4.5, 
and both they and their afterglows are easily detectable out to z > 8. Thus GRBs show great 
promise as cosmological "yardsticks" to measure the rate of expansion of the universe over time, 
and therefore the properties of dark energy (i.e., CIm, ^de, wq, and Wa). 

Measurements of four observational quantities are needed in order to use a GRB as a standard 
candle": (1) the observed peak energy E*^^^ of the burst spectrum in iyF,y (2) the energy fluence 
Se of the burst, (3) the jet break time ijet of the burst afterglow, and (4) the redshift z of the 
burst. We describe a concept for a possible MIDEX-class mission dedicated to using GRBs to 
constrain the properties of dark energy that would obtain these quantities for > 800 bursts in the 
redshift range 0.1 < z < 10 during a 2-year mission. This burst sample would enable both Qm 
and wq to be determined to ±0.07 and ±0.06 (68% CL), respectively, and Wa to be significantly 
constrained. Probing the properties of dark energy using GRBs is complementary (in the sense 
of parameter degeneracies) to other probes, such as CMB anisotropies and X-ray clusters. GRBs 
are complementary to Type la supernovae because of the broader redshift range over which GRBs 
occur, and the ease with which they and their afterglows can be observed. 

The MIDEX-class mission concept calls for three small spacecraft placed at L2 in the Earth-Sun 
system. One spacecraft, a minisatellite, would carry an array of X-ray CCD cameras very similar 
to the Soft X-Ray Camera (SXC) on HETE-2 and a set of large-area Nal detectors. The array 
of X-ray CCD cameras would have a sensitivity ~ twice that of the Wide-Field X-Ray Monitor 
(WXM) on HETE-2 . The field of view (FOV) of the array would be 27r steradians, 10 x 
greater than that of the WXM. It would localize bursts to l"accuracy. The set of Nal detectors, 
combined with the array of X-ray CCD cameras, would determine E^^^ and Se to 10% accuracy 
(90% CL) over the 1-1000 keV energy range. The other two spacecraft-both microsatellites- 
would carry identical, high throughput, short (80 cm) focal length X-ray telescopes for monitoring 
the X-ray afterglows of each of the bursts and determining tj^t for them. The mission concept 
incorporates six dedicated, ground-based 2-meter class telescopes with integral field spectrographs 
(IFS) distributed in longitude around the Earth for promptly observing the optical afterglows of 
the bursts and determining their redshifts. 
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1. Precursor Observations 

Recently, Amati et al. (2002) discovered a 
relation between the peak energy -Epeak of the 
burst spectrum in vF^ and the isotropic-equivalent 
energy £^iso of the burst in the source rest frame, 
using 10 BeppoSAX GRBs with known redshifts. 
Lamb et al. (2004; see also Sakamoto et al. 
2004, 2005) confirmed this relation for GRBs and 
found evidence that it extends down to X-ray- 
rich GRBs (XRRs) and X-Ray Flashes (XRFs) 
and therefore spans five decades in Eiso, 

Very recently, Ghirlanda ct al. (2004a) pro- 
vided further confirmation of this relation and 
discovered that an even tighter relation exists be- 
tween -Epeak and Ej, where E^ = (1— cos ^jet)-Eiso 
(see Figure 1), where ^jet is inferred from the 
jet break time tjct- Still more recently, Liang & 
Zhang (2005) have shown that an equally tight 
relation exists directly among .Epeak, Eiso, and 
tjct, which places the -Epeak — Ej relation on an 
entirely empirical footing. 

Using the -Speak — Ej relation, Ghirlanda et 
al. (2004b) and Firmani ct al. (2005) have shown 
that GRBs can be used as standard candles for 
cosmology. In, particular, they have shown that 
the uncertainty in Ejso using this method is cur- 
rently about 0.25 dex, which is only twice the 
current uncertainty in the luminosities Lsn of 
Type la supernovae (SNe) after using the anal- 
ogous relation between Lsn and the rate of de- 
cline of the SN light curves Riess et al. (2004). 
These results suggest that, although the use of 
GRBs as standard candles is less than a year old, 
while the use of Type la SNe is more than a half 
century old, GRBs hold great promise for cos- 
mology. GRBs arc also complementary to Type 
la SNe because of the broad redshift range over 
which GRBs occur, and the ease with which they 
and their afterglows can be observed. 



The "Gold" GRB sample used by Ghirlanda 
et al. (2004; see also Ghiselfini et al. 2005 
and Firmani ct al. 2005) consists of 15 GRBs 
for which -Bpeak, Eiso, ijet, and z are well deter- 
mined. Ghirlanda et al. (2004b) and Firmani et 
al. (2005) show that, using this small sample of 
GRBs, interesting constraints can be placed on 
i^M, ^DE, Wo, and wi, assuming w{z) = wq+wiz 
(see Figures 2-4). 

Precursor observations are needed to deter- 
mine E^^^ Se, ijet) and z for a much larger sam- 
ple of GRBs, and to use this much larger sample 
of bursts to confirm that (1) most GRBs satisfy 
the -Bpcak — -Eiso — ^jet relation, (2) this relation 
extends down in -Eiso and -Epcak to XRFs, (3) the 
systematic errors in this relation are small, and 
(4) X-ray afterglow observations can accurately 
determine tjet for most bursts, including XRFs, 
for which tjet can be as large as ~ 20 days. 

Confirming each of these four items will re- 
quire (1) the prompt localization of many more 
GRBs and XRFs; (2) the measurement of their 
observed peak energy 

and fiuence Se; (3) 
the determination of their jet break time tjet; and 
(4) the measurement of their redshift z. HETE- 
2 and Swift, working together, can do (1) very 
well, enabling ground-based telescopes to do (4). 
HE TE- 2 is ideally suited to do (2), whereas the 
Swift BAT (which has a relatively narrow en- 
ergy band: 15 keV < .E < 150 keV) cannot. 
Swift XRT is ideally suited to do (3), whereas 
HETE-2 cannot. However, a scientific partner- 
ship between HETE-2 and Swift can obtain the 
necessary precursor observations. 

This scientific partnership works in two ways: 

(1) a few bursts will be observed by both satel- 
lites, since HETE-2 points anti-sun and Swift is 
biasing its spacecraft pointing direction toward 
the anti-sun direction as much as it can; and 

(2) Swift can slew to the 20-25 bursts per year 
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that HETE-2 localizes and for which it measures 
-^pcak ^E- The former has already happened 
for GRBs 050215b (Sakamoto et al. 2005b) and 
050401 (Atteia 2005; Barbier et al. 2005); the 
latter has already happened for GRBs 050209 
(Kawai et al. 2005) and 050408 (Sakamoto et al. 
2005c; Wells et al. 2005). 

Of the 20-25 bursts per year that HETE-2 
localizes, ~ 7 arc GRBs, « 7 arc XRRs, and ~ 
9 are XRFs. This implies that, in three years 
of partnership with Swift, HETE-2 will localize 
and determine E^^^ and Ei^o for Ri 21 GRBs 
and XRRs, and 27 XRFs, giving a total of 
^ 70 bursts of which it can reasonably be ex- 
pected that ~ 55 will have all four quantities 
measured. A few additional bursts localized by 
Swift will have their -Epg^j^ and Se measured by 
Konus-Wind. Thus precursor observations can 
provide 60 additional "Gold" bursts - nearly 
doubling the size of the current sample of bursts 
with known redshifts, quintupling the size of the 
current sample of "Gold" bursts, and increas- 
ing the number of XRFs with known redshifts 
from two to 15-20 [XRFs 020903 (Soderberg et 
al. 2004) and 040701 (Kelson et al. 2003) are 
the only XRFs with accurately known redshifts; 
redshift constraints exist for two others: XRFs 
020427 (Amati et al. 2004) and 030723 (Fynbo et 
al. 2004)]. The full sample of 70 "Gold" bursts 
will enable a much more stringent test of the de- 
gree to which GRBs obey the Ep^nk ^ ^iso ~ ijet 
relation, particularly XRFs. In particular, con- 
firmation that XRFs obey the i^pcak ~ -E-iso — ^jet 
relation would make it possible to use them to 
calibrate the relation (see below) , and to increase 
the power of the constraints on wq and Wa (also 
see below). 

2. Error Budget 

Statistical Errors: The primary sources of 



statistical errors are the uncertainties with which 
the spectral parameters E'^^^, Se, and the jet 
break time tjct, can be determined. In addi- 
tion, there are sample variances in the -Epeak — 
Eiso and -Epeak — Ej relations of sizes 0.2 and 
0.1 dex, respectively. Here we assume Gaussian 
distributions perpendicular to both the best-fit 
-E'pcak ~ ^iso and the best-fit E'pcak ~ ^7 relations 
with widths a = 0.4 and a = 0.15, respectively. 
These widths accurately characterize the current 
statistical plus sample variances in these distri- 
butions. We assume 20% errors in both E^eak 
and tjet, and 10% errors in the fluence Se- The 
accuracies with which the parameters (intercepts 
C, slopes s, and intrinsic widths a) of the -Epcak" 
Eiso and E^f,^]^-E^ relations can be determined 
scale as I/x/IVgrb, where Nqy^b is the number 
of bursts in the sample used. The uncertainties 
in the determination of J7m, ^^de, wq, and Wa, 
therefore scale similarly. 

The scientific partnership between HETE- 
2 and Swift discussed in the previous section 
can produce significantly better determinations 
of the spectral parameters -Epeak, -Eiso 5 and E^ 
of the burst, and the jet break time ijet, espe- 
cially for the important XRFs. The mission con- 
cept that we propose would produce much bet- 
ter determinations of these quantities. Never- 
theless, here we assume the previous accuracies 
with which these quantities can be determined 
(i.e., 10% errors in Se and 20% errors in E^^ 
and tjct)- In this sense, the results that we report 
in this white paper are very conservative. 

Systematic Errors: The use of GRBs to con- 
strain the properties of dark energy is in its in- 
fancy compared with many of the other approaches 
that are being used and proposed. We acknowl- 
edge that in this situation, the risk that sys- 
tematic errors may prevent this approach from 
reaching the accuracy necessary to tightly con- 
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strain wq and especially, to remove the degener- 
acy between wq and Wa that exists if CMB priors 
arc used (as wc do) may be greater than in some 
other methods. However, there are reasons to 
believe that GRBs may suffer less from system- 
atic errors than some other approaches. 

Friedman & Bloom (2005) have discussed a 
number of possible sources of systematic errors 
in using GRBs as standard candles. Ghirlanda 
et al. (2005) and Firmani et al. (2005) have ad- 
dressed many of these. In addition, very re- 
cent work by Liang &; Zhang (2005) has shown 
that the -Epeak — Eiso — Ej relation can be put 
on an entirely empirical footing by writing it in 
terms of E'pcak — ^iso — ^jct- This approach has 
the advantage of making the relation explicitly 
model-independent, and eliminates the need to 
marginalize over the density n of the circumburst 
medium. 

Perhaps the most worrying kind of system- 
atic error for relative standard candles is one 
that would mimic the effect of dark energy; e.g., 
calibration errors that would produce a bump 
in the inferred intrinsic brightness over the red- 
shift range 0.6 < ^ < 1 (i.e., Az/z ^ 40%, 
corresponding to energy or wavelength ranges 
AE/E « AA/A « 40%). GRBs may be less 
prone to this particular kind of systematic error 
than other standard candles for several reasons. 
First, the determinations of z and tjet are not 
subject to this kind of systematic error. Second, 
the spectra of GRBs are smoothly broken power 
laws that span a broad energy range and that 
have a peak energy in uFi, that is broad. Third, 
while the detector response matrices (DRMs) of 
gamma-ray detectors (e.g, Nal) have an absorp- 
tion feature aX E 37 kev, the DRMs are not 
diagonal: because of Compton scattering and 
other effects, a beam of mono-energetic photons 
produces counts in a range of energy loss bins. 



which smooths the detector response as a func- 
tion of energy. These properties of GRB spectra 
and detectors lessens the probability that sys- 
tematic errors of the kind that would mimic the 
effect of dark energy can occur. 

Self- Calibration: The Ep^ak — -£'iso ~ ijet re- 
lation is cosmology dependent because Eiso de- 
pends on cosmology (Ghirlanda et al. 2004b). 
However, the fact that the 15 bursts in the cur- 
rent "Gold" sample of GRBs span the broad red- 
shift range z = 0.17 - 3.2 is evidence that the 
relation does not evolve with redshift; i.e., it is 
independent of z. Assuming that the £^peak — 
-E'iso — tjet relation is independent of self- 
calibration approach can be used to eliminate 
the dependence of its slope on cosmology. In this 
approach, either a sample of XRFs at low red- 
shifts or a sample of GRBs in a narrow redshift 
range at high redshifts can be used to determine 
the relation, independent of cosmology. Figure 
5 illustrates the latter approach. It shows that, 
using this method, the slope of the relation can 
be calibrated to < 1% accuracy, independent of 
cosmology. Comparisons of the parameters that 
define the relation for different redshift intervals, 
and for bursts at low redshift, can be used to 
check that the slope of the relation is indeed in- 
dependent of cosmology. Figure 6 illustrates the 
much smaller sizes of the confidence regions that 
result when the slope of the Ep^ak — Ej relation 
is known, independent of cosmology. 

Lastly, the mission concept that we are pro- 
posing would provide a data base of > 800 bursts 
with well-determined spectral parameters, jet 
break times, and redshifts that will be an order of 
magnitude larger than the sample of such bursts 
that will be produced by the scientific partner- 
ship between HETE-2 and Swift . This large 
sample of "Gold" bursts will enable a search for 
the best empirical relationship among the spec- 
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tral parameters, the jet break time of the after- 
glow, and the redshift thus tightening the con- 
straints on the properties of dark energy. It 
win also provide an unprecedented data base for 
gaining a new understanding of GRB jets - and 
therefore minimizing - possible sources of sys- 
tematic error. 

3. Constraints on Nature of Dark 
Energy 

Scaling from the burst rate seen by BATSE, 
BeppoSAji , and HETE-2 , the mission concept 
we are proposing would localize and obtain the 
spectra for k. 800 GRBs per year, or 1600 dur- 
ing a 2-year mission, of which 200 per year, or 
400 during a 2-year mission, would be XRFs. 
The percentage of HETE-2 bursts for which red- 
shifts have been obtained during the past three 
years is ~ 30%; however, most of these bursts 
had T — 10' WXM localizations, requiring rela- 
tively wide-field optical cameras to identify the 
optical afterglow, a prerequisite to determining 
the redshift of the burst from absorption lines 
in the afterglow spectrum. In contrast, the mis- 
sion concept that we are proposing would obtain 
and disseminate arcsecond burst localizations in 
real time. We therefore estimate that redshifts 
will be obtained for 50% of the bursts that the 
mission localizes, a percentage that we regard 
as quite conservative. This gives 400 bursts per 
year, or 800 bursts in a 2-year mission, with well- 
measured values of -E^peak' ^E, ijet5 and z. 

We compute the constraints that GRBs can 
place on the properties of dark energy as fol- 
lows. First, we generate a catalog of 800 GRBs 
with spectral properties matching the distribu- 
tions of the spectral properties of GRBs seen by 
Beppo^KX. and HETE-2 and a redshift distribu- 
tion matching the inferred star formation rate as 
a function of redshift (see Figure 7). We param- 



eterize the dark energy equation of state w{z) as 
a function of redshift by the phenomenological 
expression (hinder 2003) w{z) = wq + {\ — a)wa, 
where a = 1/(1 + z) is the scale factor of the uni- 
verse. In constructing the catalog of simulated 
bursts, we assume Qm = 0.30 and Ode = 0.7, 
wq = —1, and Wa = (i.e., a spatially flat, 
vanilla ACDM cosmology). 

We determine the best-fit cosmology and 
68%, 90%, and 99% confidence regions by com- 
paring the luminosity distance d^iz) for various 
cosmologies with the distribution of d^^^ values 
for the 800 GRBs in our simulated burst sample, 
using as the statistic. In doing so, we employ 
two CMB priors: + ^DE = 1 (a flat universe) 
and Q.uh'^ = 0.14. Because Planck is expected to 
determine the value of the latter to < 1%, and to 
simplify the computations, we assume that both 
priors have zero uncertainty. Given these priors, 
d^iz) is a function only of Hq, wq, and Wa- We 
consider three cases: 

1. We first assume that wq = —1 and Wa = 0. 

The only free parameter is then Hq. Figure 
8 shows the resulting posterior probability 
distributions for Qm and Qbe- 

2. We assume that Wa = 0. The free param- 
eters are then Hq and wq. Figure 9 (top 
panel) shows the resulting 68%, 90%, and 
99% confidence regions in the (i}M,wo)- 
plane. Figure 9 (bottom panel) is the same, 
except that we also employ the prior Hq = 
68 ± 8 km s"i Mpc^^ after the HST key 
project (Frccdman et al. 2001), but adopt- 
ing a best-fit value for Hq that is consis- 
tent with the values of J7m and J7de that 
we used in constructing our catalog of sim- 
ulated bursts. 

3. We make no assumptions about wq and Wa- 
We again employ the prior on Hq from the 
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HST key project, but to simplify the com- 
putations we assume that it has zero uncer- 
tainty (this is a good approximation, since 
the constraints derived on wq and Wa are 
much larger than this.) The free parame- 
ters are then wq and Wa- Figure 10 shows 
the resulting 68%, 90%, and 99% confi- 
dence regions in the (■u;0)'"^a)-plane. 

Table 1 gives the 1-a uncertainties in the dark 
energy parameters for each of the above cases, 
marginalized over the other parameters. The as- 
sumed values of the parameters lie within the 
68% confidence regions but not at their centers 
(unlike Fisher matrices) because we use a simu- 
lated catalog of GRBs. 

GRBs, like Type la SNe, are "relative" stan- 
dard candles, in the sense that their intrinsic lu- 
minosities are not determined; instead, their rel- 
ative brightnesses are compared at different red- 
shifts. Many GRBs occur and can be observed 
at high redshift. Using CMB priors then pro- 
vides an absolute calibration (Hu 2005). In this 
case, one is measuring AD/D and the most se- 
vere constraint on whether wq = —1 comes from 
GRBs observed at z < 0.3. This is opposite 
to the case for Type la SNe, since many Type 
la SNe occur at z < 0.3; then AHqD/HqD 
AHq/Hq, and the most severe constraint on whe- 
ther Wq = —1 comes from Type la SNe observed 
at redshifts z > 1 (sec Figure 11). 

Therefore, in order to exploit the power that 
comes from using CMB priors, a sufficient num- 
ber of bursts must be observed at very low red- 
shifts. If the precursor observations described 
above confirm that XRFs satisfy the -Epeak — 
-E-iso — tjct relation, XRFs can provide the pop- 
ulation of relative standard candles at z ~ 
that is needed: XRFs have redshifts z < 0.5 and 
by adjusting the design of the proposed mission 
concept, as many as 200-300 XRFs could be ob- 



tained at very low redshifts. Alternatively, the 
complementary redshift ranges of Type la SNe 
iO<z< 1.7) and classical GRBs (0.5 < z < 10) 
could be exploited to provide a combined popula- 
tion of relative standard candles that spans the 
entire redshift range < z < 10, which could 
provide more severe constraints on the proper- 
ties of dark energy than could either population 
alone. 

4. Proposed Mission Concept 

The mission concept (see Appendix A) is 
optimized to use GRBs as "standard candles" 
to determine the properties of Dark Energy and 
probe the early universe rather than increasing 
our understanding of GRBs and their precursors. 
In this sense, it is very different than previous 
GRB missions, including BeppoSAji, HETE-2, 
and Swift, although it will provide a great deal 
of information about the prompt and afterglow 
emission of a large number of bursts. The mis- 
sion concept that we describe also builds on the 
"lessons learned" from previous GRB missions, 
including BeppoSAji, HETE-2 , and Swift. In 
particular, it exploits the very recent discoveries 
that tight relations exist between Eiso and £?peak, 
and between and -Epeakj and it combines the 
capabilities of all of these missions in a way that 
maximizes the scientific return of the mission for 
the purpose of determining the properties of dark 
energy. 

As described more fully in Appendix A, the 
mission concept calls for three small spacecraft 
operating in near-proximity at the sun-earth L2 
point for two years with the following roles: (1) 
One Prompt Satellite (a "mini spacecraft," ^^200 
kg) to detect and fully characterize the prompt 
spectra of GRBs and XRFs over the 1-1000 keV 
band, as well as to establish their location to 
~arcsecond accuracy. (2) Two Afterglow Satel- 
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lites ("micro spacecraft," ~50 kg each) to mea- 
sure the burst X-ray afterglow and to determine 
the time of the jet break. In addition, as de- 
scribed in Appendix B, the mission concept calls 
for: (3) A network of six fully-automated, two 



and are small in diameter. There X-ray CCD 

detectors are even smaller than those that were 
flown on ASCA. The required sensitivity for this 
mission concept is achieved primarily by greatly 
increasing the number of detector modules, so 



meter-class ground-based telescopes equipped with as to view a much larger solid angle (~ 10 times 



integral-field spectrometers (IFS) that will be 
dedicated to establishing rcdshifts for a large frac- 
tion (~50%) of the localized GRBs and XRFs. 

5. Other Issues 

Risks and Strengths: This mission concept is 
based on a low-risk approach that relies on flight- 
proven technologies and the heritage of HETE-2. 
Redundancy is provided by a cluster of silicon X- 
ray CCD detectors and a set of Nal detectors on 
one satellite, and by two focusing X-ray telescope 
arrays, each on a separate satellite. The number 
of GRBs observed is large; the contribution of 
statistical errors to the uncertainties with which 
the dark energy parameters can be determined 
will be modest. The level at which systematic er- 
rors will dominate is not known, and represents 
a risk, but precursor observations can substan- 
tially reduce this risk. In addition, the properties 
of GRB spectra and gamma-ray detectors lessen 
the probability that systematic errors of the kind 
that would mimic the effect of dark energy can 
occur. 

Technology Readiness: There are no compli- 
cated requirements on the satellites or the oper- 
ations. The Prompt Satellite is of modest size; 
the X-ray CCD detectors are similar to those 
used in the Soft X-ray Camera (SXC) on HETE- 
2 and the focal plane of Chandra, and the Nal 
detectors are nearly identical to those that have 
flown on numerous GRB missions over the past 
30 years-most recently Fregate on HETE-2. The 
two X-Ray Afterglow Satellites are of very small 
size; the X-ray mirrors are of short focal length 



that of HETE-2 and Swiff) and by operating at 
L2, which permits a much greater observing ef- 
ficiency (~ 3 times that of HETE-2). The ded- 
icated ground segment of the mission relies on 
purchasing six commercially-available two-meter 
class telescopes and equipping them with effi- 
cient, state-of-the-art integral field spectrome- 
ters. Thus, the mission concept relies on flight- 
proven technologies and established ground-based 
telescope systems, and is thus "ready to go." 

Relationship to JDEM and LSST: The GRB 
Dark Energy mission is not a precursor to JDEM 
or LSST. Indeed, the significantly different pa- 
rameter degeneracies of the constraints on the 
parameters that characterize Dark Energy that 
will come from the GRB Dark Energy mission 
make it complementary to JDEM and LSST, as 
well as to Planck and any X-ray survey of galaxy 
clusters. 

Access to Facilities: All facilities that are 
needed for the baseline mission have been in- 
cluded as part of the mission. Ground-based op- 
tical follow-up observations using existing large- 
aperture (> 6 meter) telescopes to determine 
the redshifts of some fainter bursts (e.g., some 
XRFs) would increase the number of bursts be- 
yond that assumed in this white paper. 

Timeline: The GRB dark energy mission 
discussed in this white paper is not currently 
funded. The baseline mission fits within the re- 
sources (time, mass, volume, and cost) of a NASA 
Medium Explorer mission. It could therefore be 
proposed for next Midex AO and flown 4 years 
later (e.g., in ~ 2011). 



7 



REFERENCES Lamb, D. Q., et al. 2004, New Astronomy, 48, 

423 



Amati, L., et al. 2002, A&A, 390, 81 

Amati, L., et al. 2004, A&A, 426, 415 

Atteia, J.-L. 2005, private communication 

Barbier, L., et al. 2005, GCN Circular No. 3162 

Bloom, J. S., Frail, D. A., &: Kulkarni, S. R. 
2003, ApJ, 594, 674 



Liang, E., Sz Zhang, B. 2005, ApJ, submitted 
(astro-ph/0504404) 

Linder, E. V. 2003, PRL, 90, 091301 

Perlmutter, S., et al. 1999, ApJ, 517, 565 

Riess, A. G., et al. 2004, ApJ, 607, 665 (R04) 

Firmani, C, Ghisellini, G., Ghirlanda, G., & ^o^an-Robinson, M. 2001, APJ, 549, 745 

Avila-Reese, V. 2005, MNRAS, 360, LI Sakamoto, T., et al. 2004, ApJ, 602, 875 

Frail, D. A., et al. 2001, ApJ, 562, L55 Sakamoto, T., et al. 2005a, ApJ, in press (astro- 

Freedman, W. L., et al. 2001, ApJ, 553, 47 pli/0409128) 

„. , A o T^, TO o^^,- A T Sakamoto, T., et al. 2005b, ApJ, to be submitted 

Friedman, A. S. & Bloom, J. S. 2005, ApJ, 627, ' ' > > 

1 Sakamoto, T., et al. 2005c, GCN Circular No. 



Fynbo, J. P. U., et al. 2004, ApJ, 609, 962 



3189 

Soderberg, A. M., et al. 2004, ApJ, 606, 994 



Ghirlanda, G., Ghisellini, G., Sz Lazzati, D 

2004a, ApJ, 616, 331 Spergel, D. N. et al. 2003, APJS, 148, 175 

Ghirlanda, G., Ghisellini, G., Lazzati, D., & Fir- Wells, A. A., et al. 2005, GCN Circular No. 3191 
mani, C. 2004b, ApJ, 613, L13 

Ghirlanda G., Ghisellini, G., Lazzati, D. & Fir- 
mani, C. 2005, in Proceedings of the 
Fourth Rome Workshop on GRBs in the 
Afterglow Era, II Nuovo Cimento, in press 
(astro-ph/0504184) 

Hamuy, M. et al. 1996, AJ, 112, 2398 

Hu, W. 2005, ApJ, in press (astro-ph/0407158) 

Kawai, N., et al. 2005, GCN Circular No. 3013 

Kelson, D. D., et al. 2003, GCN Circular No. 
2627 



Lamb, D. Q., & Reichart, D. E. 2000, ApJ, 536, 
1 



This preprint was prepared with the AAS 
macros v5.2. 



Table 1. Dark energy parameter constraints from 800 GRBs.° 



Priors" 


C7(0de)' 


a{wo)' 




Wo = -1, Wa = 


-0.0076, +0.0142 






h = 0.68 ± 0.08, = 


±0.070 


±0.061 




nM = 0.3 




-0.2, ±0.32 


-1.12, ±0.63 



"All results assume CMB priors (i.e., Qm ± ^de = 1 and ^m^^ = 0.14), 
in addition to the priors listed. 
^AU confidence regions are 68.3% CL. 
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Fig. 1. — The rest frame Epeak^Eiso/ plane. Black open symbols represent the isotropic equiv- 
alent energy. Red filled symbols are the 15 GRBs for which a jet break was measured in their 
afterglow light curves. Blue symbols are upper/lower limits for E^. The four new GRBs are repre- 
sented as open green circles for £"130 and filled green symbols for E^. Also shown are two outliers 
(black squares) for either the Amati and the Ghirlanda correlation (filled squares). Stars are the 
two XRF with known redshift. The Amati correlation is also reported either fitting with the errors 
on both coordinates (long dashed line) and with the least square method (dashed line). The best 
fit Ghirlanda correlation (solid black line), giving a reduced = 1-33 and a slope ~ 0.7, is also 
shown with its uncertainty region (shaded area). (Figure from (Ghirlanda et al. 2005)). 
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Fig. 2. — Classical Hubble diagram in the form of luminosity-distance Dl vs. redshift z for 
Supernova la (open green circles: Calan/Tololo sample (Hamuy et al. 1996); open blue circles: 
Perlmutter et al. 1999) and GRBs (filled red circles: the 15 bursts in Ghirlanda et al. 2004b). In 
the top panel the SN la and GRBs are treated as standard candles (no corrections applied); for 
GRBs = 10^^ erg is assumed. In the bottom panel, we have applied the stretching-luminosity 
and the ii^^-ii^peak relations to SN la and GRBs, respectively, as explained in the text. Note that, for 
GRBs, the applied correction depends upon the assumed cosmology: here we assume the standard 
Q^M = 0.3, Ode = 0.7 cosmology. Both panels also show different D-l{z) curves, as labelled. /,From 
Ghirlanda et al. (2004b). 
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Fig. 3. — Constraints in the — r^oE (f^DE = ^^a) plane derived for our GRB sample (15 
objects, red contours); the "Gold" Supernova la sample of R04) (156 objects, blue contours, derived 
assuming a fixed value of Hq = 65 km Mpc~^, making the contours slightly different from Fig. 
8 of R04). The WMAP satellite constraints (black contours, Spergel et al. 2003) are also shown. 
The three colored ellipsoids are the confidence regions (dark green: 68%; green: 90%; light green: 
99%) for the combined fit of SN la and our GRB sample. For GRBs only, the minimum Xred ~ 1-04, 
is at Qm = 0.07, J^de = 1-2. From Ghirlanda et al. (2004b). 
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-2.0 -1.5 -1.0 -0.5 0.0 

Wo 

Fig. 4. — Constraints on the f^O) w' parameters entering the equation of state p = {wq + w'z)pc^, 
where p is the dark energy density, wq = —1 and w' = correspond to the cosmological constant 
r^DE- We assume a flat geometry and Qm = 0.27 (see also R04). Blue contours: constraints from 
type la SN (R04). Red contours: constraints from our GRBs, Colored regions: combined constraints 
(dark green, green and light green for the 68%, 90% and 99% confidence levels, respectively). From 
Ghirlanda et al. (2004b). 
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Fig. 5. — Illustration of the self-calibration method for the -Bpcak ~ -E'iso relation. For different 
redshift "slices" Az, centered on redshifts z = 0.7, 1.0, 2.0, and 3.0, the following quantities are 
shown: Top panel: Maximum variation 6m in the slope m of the -Epeak ~ correlation, obtained 
from surveying large ranges in the cosmological parameters Qm and JIde- The numbers next to 
each point are the number of bursts in that particular redshift slice. Bottom panel: Maximum 
variation in obtained from surveying large ranges in the cosmological parameters Q,m and Q^e- 
Also shown is the 68% confidence level. A value of Ax^ below this level ensures that the -Epeak — -E'iso 
relation found using that redshift range (i.e., that size of Az) is cosmology independent. 
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0.0 0.5 1.0 1.5 2.0 




Fig. 6. — Example of how GRBs can contribute to the determination of the cosmological param- 
eters once the — E'peak correlation can be determined in a cosmology-independent way [i.e., 
accumulating sufficient bursts at small redshifts or in a small redshift range (e.g. around z = 2)]. 
We show the contours in both the {0,m, ^^de)- plane (Ode = in main figure) and in the wq — w' 
plane (insert, a flat cosmology with JIm = 0.27 is assumed). Lines and colors are as in Fig. 2 and 
Fig. 3. (Figure from Ghirlanda et al. (2004b)). 
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Fig. 7. — Distribution of GRBs as a function of redshift in our simulated catalog of 800 bursts. The 
distribution of GRBs is proportional to the phenomenological expression for the star formation rate 
(SFR) suggested by Rowan-Robinson (2001), adopting the parameter values P = 1.2 and Q = 5.4 
(which provides a good fit to existing data on the SFR as a function of redshift) . With this choice 
of the SFR and the nominal mission design, the catalog contains 35 GRBs at 2; < 1. However, 
XRFs have redshifts z < 0.5, and by adjusting the design of the proposed mission concept, the 
catalog could contain as many as 200-300 XRFs at very low redshifts. 
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Fig. 8. — Constraints on ^Im and Ode for 800 GRBs, using CMB priors (the only free parameter 
is then Hq). Top paneh posterior probabihty distribution for i^M- Bottom panel: posterior 
probability distribution for r^oE- The dotted lines show the 68% credible regions; the dashed lines 
show the 90% credible regions. 
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Fig. 9. — Top panel: ^Im versus wq for 800 GRBs and using CMB priors. Bottom panel: the same, 
except that we also use the prior Hq = 68±8 km s~^ Mpc~^ after the HST key project (Preedman 
et al. 2001). 
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Fig. 10. — Constraints on wq and Wa for 800 GRBs, using CMB priors and the prior Hq = 68 km 
Mpc~^ after the HST key project (Freedman et al. 2001). The last prior is a good approximation, 
since the constraints derived on wq and Wa are much larger than this. However, using the correct 
prior, Hq = 68 it 8 km s~^ Mpc~^ [after the HST key project (Freedman et al. 2001)], would widen 
the contours in the direction perpendicular to their long axes. 
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Fig. 11. — Sensitivity of determination of wq, (assuming Wa = 0) as a function of redshift for 
various quantities. Both GRBs and Type la SNe are "relative standard candles," in the sense that 
their intrinsic luminosities are not determined; instead, their relative brightnesses are compared at 
different redshifts. Many GRBs occur and can be observed at high redshift. Using CMB priors 
then provides an absolute calibration (Hu 2005). In this case, one is measuring AD/D and the 
strongest constraints on whether wq = —1 come from GRBs observed at z < 0.3. This is opposite 
to the case for Type la SNe, since many Type la SNe have been observed at z < 0.3; then 
AHqD/HqD — > AHq/Hq, and the strongest constraints on whether wq = —1 come from Type la 
SNe observed at redshifts z > 1. From Hu (2005). 
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Appendices: Instrumentation Concept for 
GRB Dark Energy Investigation 
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Introduction 

As outlined in the text, for each of the 800 
GRBs in the sample, we need to measure the 
spectral peak, -Epeaki the isotropic energy flu- 
ence, -Eiso; the jet break time, tjet; and the red- 
shift, z. In the concept we describe below, sensi- 
tive silicon and sodium iodide detectors mounted 
on a "Prompt Satellite" will establish -Bpcak a-nd 
-Eiso, as well as burst location. A dedicated op- 
tical ground network will measure the redshift, 
as described in Appendix B. Compact X-ray 
Afterglow Telescopes (XATs), flying on two "Af- 
terglow Satellites" will measure tjet- Figure A.l 
depicts the interaction between the assembled 
assets for the mission concept. 

The key to measuring tjct and z is rapid 
and precise determination of astrometric posi- 
tions, allowing the optical ground network and 
XATs to acquire the afterglow while it is still 
bright. Soft X-ray Camera (SXC) modules with 
5 ;um pixels on the Prompt Satellite designed 
for precision astrometry should be capable of 
locating bursts to an accuracy of l"-3" in real 
time. This is sufficiently accurate for immediate 
ground-based IPS spectroscopy and for pointing 
a satellite-borne XAT with a narrow field of view. 

The HETE-2 SXC has selected a much higher 
yield of bright afterglows than other GRB locat- 
ing instruments. In addition, the distribution of 
redshifts in this sample is very flat out to ^; > 3 
(Berger et al. 2005). This is exactly the kind of 
sample that is required for a dark energy inves- 
tigation. 

Relative to the HETE-2 SXC, the instru- 
ment concept we propose will cover a much greater 
solid angle with much higher operational effi- 
ciency. This broad coverage is the only way 
to acquire large numbers of bright bursts. The 
HETE-2 SXC localizes ~ 12 bright bursts/year 
in a ~ 0.6 sr field of view with an operational effi- 



ciency of ~ 20%. An array of HETE-2 SXCs with 
hemispherical coverage and 80% operational ef- 
ficiency would therefore localize ~ 500 bright 
bursts/year. This is a very conservative lower 
limit on the burst rate for this instrument con- 
cept, which will have ~ 30 times the detector 
area per solid angle and will therefore be far more 
sensitive than HETE-2. If the burst rate scales 
as A^-/^ (Doty 2004), this would imply that we 
might anticipate ~ 1600 events per year. We 
adopt a conservative value of 800 per year. 

The Prompt Satellite will also carry a large 
array of sodium iodide (Nal) scintillators to ac- 
curately measure burst spectra and fluences in 
the 5-1000 keV range. The technology for such 
detectors is well-established, based on their use 
in GRB satellites from the 1960s until the present- 
most recently in the Pregate instrument on HETE- 
2. Compared to Pregate, the proposed Nal array 
will provide 75 times larger area, 3 times greater 
solid angle, and 3 times greater observing time 
per year. Based on the Pregate rate of ~ 90 
bursts/year, we conservatively estimate that the 
Prompt Satellite wiU detect ~ 1000 bright GRBs 
and XRPs per year-from which we can extract 
high S/N spectra for accurately measuring a, 
P, E'pcakj and the bolometric fluence over the 
5-1000 keV range. The SXC silicon detectors 
will provide complementary spectra in the 1-10 
keV range with the required sensitivity. In the 
following two appendices, we provide additional 
details on the space and ground segments of the 
mission concept, respectively. 
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Fig. A.l. — Schematic diagram of spacecraft locations and operation. 
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A. Instrument Concept for Space-based 
Segment— GRB Satellites at L2 

A.l. Prompt Satellite 

The Prompt Satellite is depicted in Figure 
A. 2. Prom its location at L2, the detectors on 

its seven facets view the sky continuously. Each 
facet contains 14 wide-angle (~ 140° FOV) Nal 
detectors (1750 cm^ area per facet), covering the 
5-1000 keV energy range. Because the fields- 
of-view of the Nal detectors on adjacent facets 
overlap, it is possible to establish crude (~ 6°) 
localizations from detector counting rate ratios 
(as was done for the BATSE mission). Addi- 
tional details on the Nal detectors are given in 
Table A.l. 

To refine these crude localizations and to 
extend the X-ray energy range down to ~1 keV, 
each facet also has eight large area Soft X-ray 
Cameras (SXCs) with X-ray CCD detectors (Ta- 
ble A.l). The eight detectors are equipped with 
different coded masks: 4 have one-dimensional 
masks (2 oriented in the X direction and 2 in 
the Y direction), 2 have two-dimensional masks, 
and 2 have no masks. The ensemble of eight 
SXC detectors works in concert to provide excel- 
lent burst detection sensitivities in the 1-15 keV 
band, as well as localization accuracy of 1.3" in 
radius (90% C.L.). 

The seven facets on the Prompt Satellite 
host 56 SXC cameras in total. Each of the 56 
SXC cameras contains a wafer-scale (~90 cm^) 
CCD consisting of 16x16 panes with 1152 xll52 
pixels each. The multipane approach has many 
advantages. Each charge packet moves only a 
small distance (1.2 cm at most) to the readout, 
so that charge loss from radiation induced traps 
is minimized. Bad panes, caused by manufactur- 
ing defects or by niicronicteoritcs that penetrate 
the camera entrance window, may be isolated 



and ignored: coded aperture imaging is insensi- 
tive to small gaps in focal plane coverage. Thus, 
a high manufacturing yield of acceptable wafers 
as well as a long detector life in orbit is assured. 
Each pane includes charge summation structures 
to make pixel summation fast and efficient. For 
the ID modules, 576 rows are summed to pro- 
duce a tall (2.9 mm) narrow (5 ^va) effective 
pixel shape. For the 2D and open modules, the 
sums will be 24x24, for a 120 //m square effective 
pixel. This reduces the readout load from 19 bil- 
lion pixels to 33 million. A 200 Mpix/s readout 
system will allow 4 Hz readout of the entire array 
of modules with allowance for overhead. Thus, 
the CCD detector time resolution is improved 
to 250 ms, 5 times better than for the HETE-2 
SXCs. 

A large array of SXC modules requires a 
high performance data acquisition and process- 
ing system. We have previously investigated the 
design of such a system for a possible mission of 
opportunity. It appears practical, with technol- 
ogy available in 2005, to acquire and process X- 
ray CCD data for ~ 1/ 60 the power needed with 
HETE (1992) technology. HETE-2 only provides 
~ low of power to support the SXC, but the 
larger spacecraft we envision should be able to 
provide at least lOOW to the SXC array, so we 
anticipate a data acquisition and processing sys- 
tem with at least 600 times the HETE-2 capacity 
supporting ~ 200 times the silicon area. The ex- 
tra capacity will provide better time resolution 
and better real time imaging. Additional details 
on the SXC camera and X-ray CCD detectors 
are given in Table A.3. 

The sensitivity of the Prompt Satellite for 
establishing E^^ak is depicted in Figure A.3. Both 
the scintillation detections ("DE NaT") and SXC 
cameras ("DE Si") are more than an order of 
magnitude more sensitive than flown on previ- 
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Table A.l: Scintillation Detectors for Prompt Satellite 



Energy Range 

Field of View 

Number of Modules 

Scintillator 

Readout System 

Sensitive Area 

Source Localization 

Flux Sensitivity {Wa, 8-1000 keV) 

Time Resolution 

Bright Burst Rate 

Flight Heritage 

Total Power 

Total Mass 



5 - 1000 keV 
7 sr total 
14 X 7 = 98 
Nal(Tl) 

Photomultipliers (PMTs) 

12,250 cm^ (14 X 7 = 98 units @ 125 cm^) 

6° (by projected area method) 

10-8 ergs cm-2 s"^ (3x better than BATSE & Fregate) 
10 fis 

~ 1000 per year (with measured -Epeak) 
Fregate Instrument on HETE-2 
25 W 

95 kg (50 kg Scintillators +40 kg PMTs +5 kg electronics) 



ous missions. 

Despite the large areas of the Nal and SXC 
instrument arrays, they are compact, low in mass, 
and require little power. As a result, the space- 
craft power and mass resources are modest. We 
estimate that the Prompt Satellite (a "minisatel- 
lite") will require a net of < 300 W and < 400 
kg, even allowing for ample reserves of 50% and 
100%, respectively (Table A.3). 

A. 2. X-ray Afterglow Satellite 

In Figure A. 4, we show a concept for a mi- 
crosatellite optimized to observe the X-ray af- 
terglows. In Table A. 4, we have compiled the 
power and mass resources needed for each of two 
such microsatellites. Even with a 100% reserve, 
each Afterglow Satellite requires only 60% of the 
mass of HETE-2. The power requirement per 
microsatellite is approximately 60% of the power 
requirement of HETE-2. 

Each X-ray Afterglow Satellite will carry 9 



BATSE 


1 


HETE Fregate- l'^, _ 




■ DENal""""^^ 




HEtE WXM ^ 




DE Si 





1 10 100 1000 10000 

Ep*=[keV] 

Fig. A.3. — The sensitivity curves for the Dark 
Energy Si and Nal detectors (colored lines) ver- 
sus the observed GRB peak energy E'peak- Also 
plotted are the sensitivities for BATSE (dotted 
black line) and HETE-2 (dashed black and dot- 
dashed gray lines). 
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Table A.2: Silicon CCD Detectors for Prompt Satellite 



Energy range 
Field of View 
Module Configuration 
(on 1 of 7 satellite faces) 
CCD Detectors 



Pixel Summation 
Area per CCD 
Number of CCDs 
Total CCD Area 
CCD Temperature 
Coded Mask Size 
Mask-to-Detector Separation 
Camera Resolution 

Source Localization Accuracy (90% C.L.) 
Flux Sensitivity {10a, 2-10 keV) 

Time Resolution 

Bright Burst Localization Rate 

Flight Heritage 

Total Power 

Total Mmss 



1-15 keV 

7 sr total (centered on anti-sun) 
ID cameras: 2 X, 2 Y 

2D cameras: 2-D, 2 Open 

Deep depletion silicon; 72K x 72K, with 5 /xm pixels 
(Wafer scale device, segmented as 16 x 16 = 256 

independent "panes", each with 1152 X 1152 pixels) 
576 X 1 (ID) or 24 x 24 (2D) 
90 cm2 

8 CCDs/Face x 7 Faces = 56 CCDs 
56 @ 90 cm2 = 5040 cm^ 

Passive cooling to -lOOC 
10 cm X 10 cm 
20 cm 

10. 3" for Nyquist sampling of mask pattern 
1.3"radius for 5cr detection in one XY pair (2 CCDs) 
1 s: 4.5 X 10~^ ergs cm~^ s~^ 
100 s: 4.5 X 10~^° ergs cm~^ 
250 ms 

~ 800 per year 

X-ray CCDs on ASCA, Chandra, HETE-2, k Astro-E 
100 W 

■55 kg (10 kg cameras + 15 kg (Jcctrouics) 



(3 X 3) clustered, narrow field of view X-ray Af- 
terglow Telescopes (XAT), each with its own CCD 
detector. The telescope and detector design pa- 
rameters and performance statistics for a typical 
observation (field angle of 0.5') are displayed in 
Table A.2. The XAT design is an updated ver- 
sion of the array of small telescopes flown and 
proven aboard the SAS-3 satellite in the mid- 
1970s (Hearn et al. 1976). The reflecting sur- 
face of each telescope is a paraboloid of rotation, 
coated with Ni. A single bounce at grazing inci- 
dence focuses 0.5 keV X-rays onto the CCD with 
excellent efficiency. (The efficiency at 2 keV de- 



clines by only a factor of 2 compared to its value 
near unity at 1 keV.) Each telescope is baffled to 
prevent X-ray's from field angles > 1' from reach- 
ing the detector. The 9 telescope design is light 
(< 1 kg per telescope), and the small diameter 
(8 cm), thick- walled (3mm) unnested reflectors 
should be inexpensive to fabricate. 

Each X-ray Afterglow Satellite will reach a 
limiting sensitivity of 2 x 10^^ erg cm~^ in 4 
ksec. We take this flux limit as a conservative es- 
timate for the required sensitivity to detect late 
GRB jet breaJi; times. A flux ~ lOx greater was 
measured by the Swift XRT for GRB 050408 at 
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Table A. 3: Resources for GRB Prompt Satellite 



Item 


r^ower [yy ) 


iViass l^Kgj 












io ^ rjiccxromcs 


OClliLllidjLlUli l^ClcCLUl o 




l^n — NaT 


(^iz,zoU cm ) 




4U — J^iVl ± S 






5 — Electronics 


/itTjlTjUQc l^OIlXrOl oybtt:!!! 


10 




I iVlUlilCiilj Lllil W llccio , kJ Udi ti ctL^JvtJi G J 






-L ill LlO tC-i iJlCL/ ti '^iilL/O I 








5 


5 


(Electronics) 






Power System (Batteries, Converters 


5 


10 


Solar panels) 






X-Band System 


20 


1 


VHF System 


4 


1 


Propellant (Cold gas) 




5 


Sums 


174 


182 


Margin 


87 (50%) 


182 (100%) 


TOTAL (w/ Margin) 


261 W 


364 kg 



its jet break (t ~ 5 days) after > 10 ksec of in- 
tegration. Figure A. 5 shows the sensitivity of an 
XAT (solid line) relative to the sensitivity of the 
Swift XRT (dotted line). For relatively bright 
sources (fluence > 2x 10~^^ erg cm~^ s~-^), these 
small XAT instruments achieve the same signal- 
to-noise determination as docs the Swift XRT, 
yet the required exposure is 2.5 x shorter. For 
fainter sources, or for sources in crowded fields, 
the required exposure is 5 — 8x shorter than for 
the Swift XRT due to the tighter PSF of the 
XAT (10.3" HEW versus 18" PSF^). 



http: / /www. swift .psu.edu/xrt /details 



A. 3. Launch Vehicle Requirements 

As shown in Figure A. 6, the Prompt Satel- 
lite and the two Afterglow Satellites are easily ac- 
commodated within the fairing of a Delta 7425, 
the same rocket used to launch the WMAP mis- 
sion to L2 (Bennett et al. 2003). Based on the 
mass estimates for the three satellites (Tables 
A. 3 and A. 4), the Delta 7425 launch mass ca- 
pability to L2 of 800 kg provides ample (>100% 
mass margin) for the MIDEX mission we envi- 
sion. 

A. 4. Data Rate 

The data handling strategy for our mission 
combines concepts from RXTE and HETE-2. An 
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Fig. A. 2. — Schematic diagram of the GRB Prompt Satelhte. 
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RXTE-style bit code will represent raw event 
data from the scintillation detectors, with 32 chan- 
nels of energy resolution and 20 ms time resolu- 
tion. This code averages 2 bits/event. For the 
silicon detectors, the events are very sparse in 
the spatial/spectral/temporal resolution space. 
A run length encoding with 20 bits / event should 
be nearly optimal for these. The silicon detec- 
tors detect about an order of magnitude fewer 
photons than the scintillation detectors, so the 
total bit rates for the two detection systems are 
approximately equal. 

Nearly all of the time, the diffuse X-ray back- 
ground will dominate the event rate. A very in- 
tense burst may dominate the event rate for a few 
seconds, and a solar radiation storm may blind 
the detectors for awhile, but the diffuse X-ray 
rate determines the average raw data rate over 



time. For the encodings above, the average data 
rate is about 500 kbit/s. A 128 Gbit data store 
will therefore be able to hold ~ 3 days of data on 
board to support ground queries for more detail. 

Because 500 kbit/s is much larger than our 
raw downlink rate of 30 kbit/s, we will need to 
be selective in what we send to Earth. We will 
bin the raw data from each face of the local- 
ization/spectra spacecraft into 16 channels with 
1/4 second time resolution. With light data com- 
pression, this stream will be 3.6 kbit/s. This will 
constitute a standard survey data set. For most 
bursts, this should be sufficient to determine flu- 
ence and spectral parameters. 

The imaging software on board will conduct 
a spatially resolved search for bursts: we will 
devote 2 kbit/s to the results of this search. The 
on board burst trigger algorithm will provide 1 
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Fig. A. 4. — Schematic diagram of Afterglow Satellite (1 of 2) with its X-ray Afterglow Telescope 
(XAT) array. 



kbit/s of diagnostic information. Each of the 
two Afterglow Satellites will send 600 bit/s to 
the Prompt Satellite for relay to Earth. 

The above, plus spacecraft housekeeping and 
data overhead, will use 10 kbit/s, or 1/3 of the 
available downlink. The other 2/3 will be used 
to respond to ground queries for more detail. In 
most cases, we expect that all of the event data 
from a burst will be requested. For the bright- 
est/longest bursts, an RXTE-like binning and 
recoding engine will reduce the data mass to a 
manageable amount. 

For the afterglow satellites, 10 photon event/s 
will yield 320 bit/s. The remaining 280 bit/s will 
be housekeeping and overhead. In most cases the 
afterglow intensity will be much less than 1 pho- 



ton/s at the jet break time so this should provide 
adequate dynamic range. 



A. 5. Navigation and Station-Keeping 

For simplicity, the Delta 7425 payload bus 
will transport the "mother" mini spacecraft (aka 
Prompt Satellite) and the accompanying two 
"daughter" micro spacecraft (aka Afterglow Satel- 
lites) to the vicinity of the Earth-Sun L2 La- 
grangian point where they will be injected into a 
so-called halo (or Lissajous) orbit of co. 100,000- 
200,000 km radius. Such orbits (Farquhar 1970) 
lie close to a plane oriented at ~72° from the 
ecliptic. Once a halo orbit has been achieved, 
the daughter spacecraft will separate at low rel- 
ative velocities so as to execute Lissajous orbits 
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Fig. A. 6. — Schematic diagram of the GRB Prompt SateUite and 2 GRB Afterglow Satelhtes in 
the Delta 7425 rocket faring. There is ample clearance for the 3 small satellites to fit within the 
faring. Since the Delta 7425 can launch 800 kg to L2, it can easily launch the combined payload 
mass (see Tables A. 3 and A. 4), even with a mass reserve of more than 100%. 
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Table A.4: Resources for Afterglow Satellites (per Satellite) 



Item 


x^ower yy\ ) 


iviass l^Kg^ 


L/LD Detectors 


1 


2 


ieiescope Array 




10 


Attitude Control System 


1 K 


1 n 
iU 


(Momentum wheels, Star 






rracKers, ±nrusters, H/iectromcs j 






OUlllllldliU o6 J-Zata XldilLlliil^ 


o 


K 

o 


(Electronics) 






Power System 


4 


8 


(Batteries, Converters, 






Solar panels) 






VHP System 


4 


1 


Propellant (Cold gas) 




2 


Sums 


29 


38 


Margin 


14.5 (50%) 


38 (100%) 


TOTAL (w/ Margins) 


44 W 


76 kg 



(Barden &; Howell 1998) relative to the mother 
spacecraft. None of the spacecraft will be in sta- 
ble orbits — they will separate from each other, 
and from the vicinity of the L2 point, along an 
"unstable manifold" of the halo orbit at an expo- 
nentially increasing rate unless appropriate ve- 
locity corrections are made. The mother space- 
craft will contain an X-band transceiver, syn- 
chronous transponder, and high-gain antenna, 
allowing it to be tracked directly from Earth — 
utilizing round-trip time delay measurements to 
determine range, and two-way Doppler tracking 
to measure line-of-sight velocity, right ascension, 
and declination. Daily tracking as part of sci- 
ence data acquisition contacts at ~100 m accu- 
racy will be sufficient to permit velocity correc- 
tions to be computed and executed to maintain 
the halo orbit with low fuel consumption and 
slow drift orthogonal to the orbit plane (Col- 
lange & Leitner 2004). Since the total delta-V 



requirements will be < 5 m s"^ yr~^ (Parquhar 
Sz Dunham 1990), a small cold gas system (< 5 
kg propellant for the mother-sat, <2 kg propel- 
lant for each of the daughter-sats) will suffice 
for the nominal 2 year mission. We propose the 
use of a simple low-gain phased-array antenna 
on the daughter spacecraft to determine the di- 
rection to the mother spacecraft with the VHP 
uplink; 4-6 whip antennas at different places on 
the spacecraft body will feed a coherent multi- 
channel digital receiver. The phase difference 
between the antennas will indicate the direction 
of the mother spacecraft. (With antennae sepa- 
rated by less than a wavelength this is more like 
classic 1940's radio direction finding than a 21st 
century phased array, but it should suffice for our 
requirements.) This knowledge may also be used 
to phase the antennas for maximum satellite-to- 
satellite link margin. The daughter spacecraft 
will relay their data through the mother space- 
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Table A.5: X-ray Afterglow Telescope (XAT) De- 
sign Parameters 



Table A. 6: X-band Downlink Budget 



Focal Length (F) 
Mirror Length 
Mirror Outer Diam. 
Mirror Thickness 
Number of Telescopes 
per Satellite 
Mass per Telescope 
Telescope PSF 
Effective Area 
Detection Element 
Pixel Scale 

Field of View (radius) 
Sensitivity (in 4 ksec) 



80cm 
60cm 

8cm (F/10) 
3mm 

9 

0.9 kg 

10.3"HEW, 0.5 keV 
286 cm2, 0.5 keV 
5x5 micron pixel 
1.3'7pixel 

60 arc second 
2x10"^"^ erg/cm/s 



Data type 



kbit/s 



Afterglow data relay 1.2 

GRB Survey 3.6 

Soft imaging search 2.0 

Trigger status 1.0 

Housekeeping 1 . 

Overhead 1.2 

Detailed analysis 20.0 



TOTAL 



30.0 



craft, and cannot therefore be reliably tracked 
from the Earth. Instead, their range, range rate, 
and bearing from the mother spacecraft will be 
determined using their VHF (2-meter, 140 MHz) 
transceivers. Each daughter will possess a single 
dipole antenna, a VHF transceiver, and a syn- 
chronous transponder, whereas the mother will 
employ a 3 dimensional antenna array, e.g., sim- 
ple 1/4- wave VHF dipoles at each apex of an 
equilateral pyramid, each with its own receiver. 
The physical antenna can be spring-deployed (tape 
measure) elements of the type used on the HETE 
satellite, for which the VHF antenna was liter- 
ally a ~ 50 cm long spring tape measure that 
uncoiled in orbit. The mother spacecraft will 
transmit on a single dipole and its receiver out- 
puts will be heterodyned to baseband and cross- 
correllated to extract the relative phase delay 
(±10^'^ rad) of each baseline. Since the antenna 
spacing is comparable to the carrier wavelength, 
aliasing will not occur and the direction from the 
mother to each daughter will be determined to 



~1 milliradian accuracy. Combined with tradi- 
tional VHF ranging and range-rate extraction, 
the relative location of each spacecraft can be 
estimated, and velocity corrections can be com- 
puted to maintain each daughter in the vicinity 
of the mother. 



A. 6. Satellite Telecommunication 

The constellation of GRB Prompt ( "mother" ) 
and two Afterglow ("daughter") Satellites will 
circle about L2 in a halo orbit, of radius 100,000- 
200,000 km. All ground communication will pass 
through the Prompt Satellite over an X-band (10 
GHz) link, utilizing two small (~ 2000 cm^), flat 
phased array antennae on the sun- (and earth-) 
facing surface. The phased array antennae will 
permit continuous communication to earth with- 
out requiring a steerable dish or a re-orientation 
of the Prompt Satellite, which would disrupt GRB 
science data collection. The beam of the phased 
array antennae will be sufficiently narrow (< 2°) 
to avoid interference from solar RFI, yet broad 
enough to illuminate the entire earth (viewed 
from the L2 halo orbit, the earth subtends ~ 0.5° 
and is offset ~ 4° from the sun). In addition, the 
Prompt Satellite will communicate with the Af- 
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Table A. 7: X-band Telemetry Link Margins for a 5m Diameter Ground-Based Antenna 



X-band Prompt Satellite-to-Ground 



Eb/N at le-6: 
Modulation: 
Spectrum utilization: 
Spacecraft TX power: 
Distance: 
Bit rate: 
Link margin: 



7.6dB 
BPSK+l/2Viterbi 
0.5 
5W 

1,504,0000 km 
30 kbps 
4.13 dB 



X-band Ground-to-Prompt Satellite 



Eb/No at le-6: 
Modulation: 
Spectrum utilization: 
Ground station TX power: 
Bit rate: 
Link margin: 



14.2 
CPFSK 
1 

50W 
113,600 bps 
5.20 dB 



One XAT Microsatellite 
Two XAT Microsatellites 
Swift XRT 
Chandra AGIS-S 




10' 10" 10' 
Integration Time [seconds] 



Fig. A. 5. — The XAT sensitivity (colored lines) 
versus exposure time. The red line refers to 
a single XAT microsatellite viewing a given 
source, while the blue line refers to both XAT 
microsatellites simultaneously viewing the same 
source. Also plotted is the sensitivity for the 
Swift XRT (dotted hne) and for Chandra ACIS- 
S (heavy black line). 



terglow Satellites via 2-way VHF link. 

Table A. 8: VHF Telemetry Link Margins 
VHF Satellite-to-Satellite 



Frequency: 
Spacecraft TX power: 
Bit rate: 
Link margin 
(@ satellite-to-satellite 
distance) 



137.96 MHz 
1 W 
600 bps 
27 dB @ 10^ km 
7 dB @ 10^ km 



The VHF transceiver system described in 
"Navigation and Station-Keeping" will also serve 
to transmit the GRB burst alerts and ~ arcsec- 
ond localizations from the Prompt Satellite to 
the Afterglow Satellites. (An automated sched- 
uler on the Prompt Satellite will establish which 
of the two Afterglow satellites will be assigned 
to a given alert.) In addition, the VHF link will 
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be utilized to transfer the X-ray afterglow obser- 
vation data from the Afterglow Satellites back 
to the Prompt Satellite, for re-transmission to 
earth. We estimate that the satellite-to-satellite 
separations can be easily maintained to <1000 
km, even if the station-keeping thruster events 
are only scheduled a couple of times per month. 
To assess the link requirements for this VHF sys- 
tem, we have scaled from the measured proper- 
ties of the HETE VHF burst alert system. (NB: 
Although HETE is is low earth orbit and com- 
municates to ground receivers, the link distances 
are typically 1000-2000 km, which is comparable 
to the satellitc-to-satellite distances near L2 that 
wc arc considering. Furthermore, the VHF RFI 
environment is much more favorable at L2, since 
VHF RFI is dominantly terrestrial in origin. ) 

B. Instrument Concept for 
Ground-Based Segment — Dedicated 
Telescope Network 

B.l. Integral Field 
Spectrometer-equipped Telescopes 

A 2.4-meter Ritchey- Chretien telescope (such 
as production units from EOST^) will provide a 
small 30 arcsecond diameter field of view to an 
lenslet array coupled to optical fibers. An off- 
the-shelf 2.4m telescope should be able to slew 
and set to 3 arcsecond accuracy. Focus and point- 
ing will be monitored by pointing to a table of 
known standard stars during clear night opera- 
tions. 

The lenslet array will be a set of hexagonal 

elements, each subtending 0.5 square arcseconds 
(approx 0.7" on a side). We require 800 lenslets 
to cover a square 20"x20" spectroscopic FOV at 



^http: / / www.eostech.com/two_point_four_meter.php 



Cassegrain focus. The output of the fibers will 
feed a double-spectrograph, similar in design to 
the SDSS grism spectrographs. The fibers will 
be aligned and collimated to provide an average 
of 5 pixel separation between individual spectra, 
and 2000 pixels on each blue and red side (4000 
pixels total in each spectrum, with minimal over- 
lap near the dichroic). 



Table B.l: Sensitivity Estimate for Spectr- 
graphic Redshift Limits for Ground-Based IFS- 
equipped Telescopes 



ruAB 


Bright Sky S/N 
19 mag/arcsec^ 


Dark Sky S/N 
21 mag/arcsec^ 


16.0 


75.3 


78.9 


17.0 


43.4 


48.2 


18.0 


23.0 


28.2 


19.0 


11.0 


15.4 


20.0 


4.8 


7.6 


21.0 


2.0 


3.4 



If we assume pixels which are approximately 
100 km/s in dispersion, and an effective aper- 
ture of 1.0 m^ (25% end-to-end efficiency) then 
a single 900 second exposure will yield approxi- 
mately the S/N values given in Table B.l for two 
extremes in sky conditions. 

The spectrographs will have complete wave- 
length coverage from 3100-9600 Angstroms. This 
covers the absorption line doublet of Mgll over 
the full redshift range of 0.11 < z < 2.4. In 
addition, it will cover the CIV doublet absorp- 
tion from 1.0 < z < 5.0, and Damped Lyman-a 
absorption from 1.6 < z < 6.5. 

To detect Mgll and CIV robustly, one would 
require an approximate S/N of 5 per pixel. So in 
dark-sky and bright-sky (3/4 moon) conditions, 
a single 900 second exposure is sufficient down 
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Table B.2: Budget Breakdown for Dedicated 
Support Telescope k, IFS Instrument Configu- 
ration 

$2.0M 2.4 meter automated telescope: 

built, tested and delivered 
$1.5M Optical double-sided spectrograph 

with fiber /lenslet array. 
$0.3M Estimate for site lease and minimal 

observatory support over 2 years. 
$1.0M Data reduction, z determination 

and quick response support 

over the 2 year mission. 
$4.8M TOTAL per site 



to ruAB = 20th magnitude. To accurately detect 
these Mgll and CIV redshifts down to m^B=21 
will require 90 minutes of exposure in bright time 
and 30 minutes in dark lunar conditions (see Ta- 
ble B.l). The Damped Lyman- alpha feature is 
much easier to detect, and should be recovered 
in any spectra with S/N > 2.0 per pixel. 

To have complete sky coverage, we would re- 
quire 5 or 6 identical telescopes adequately sepa- 
rated in longitude. Candidate sites include: Ca- 
nary Islands, Chile, Southwestern USA, Hawaii, 
Australia, and Southern Africa. Each of the 
telescopes would be built at an established as- 
tronomical site with good infrastructure. Al- 
though the telescopes will be designed to run 
autonomously, lost observing time could be min- 
imized if observatory staff could respond in a rea- 
sonable time if trouble arose. In Table B.2, we 
give the estimated costs per dedicated observa- 
tory site for a two year mission. Since the total is 
$4.8M per site, establishing 6 sites would require 
~$30M, which would comprise ~10% of the total 
budget for the MIDEX mission. 
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